Factorization and NLO QCD correction in e+e — > J/ip(ip(2S)) + Xco at B Factories 
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In nonrelativistic QCD (NRQCD), we study e + e~ — > J/ip(ip(2S)) + XcO at B factories, where 
the P-wave state XcO is associated with an S-wave state J/ip or ip(2S). In contrast to the failure 
of factorization in most cases involving P-wave states, e.g. in B decays, we find that factorization 
holds in this process at next to leading order (NLO) in a s and leading order (LO) in v, where 
the associated S-wave state plays a crucial rule in canceling the infrared (IR) divergences. We also 
give some general analyses for factorization in various double charmonium production. The NLO 
corrections in e + e~ — > J /tj)[ip(2S)) + \c0 at ^fs — 10.6 GeV are found to substantially enhance the 
cross sections by a factor of about 2.8; hence crucially reduce the large discrepancy between theory 
and experiment. With m c = 1.5GeV and fi — 2m c , the NLO cross sections are estimated to be 
17.9(11.3) fb for e + e~ — > J/tp(i^(2S)) + Xco, which reach the lower bounds of experiment. 

PACS numbers: 13.66.Bc, 12.38.Bx, 14.40.Gx 



The production of double charmonium in e + e an- 
nihilation at B factories 0, S, i, S i, i, S @, @, E 



I 111 . Il2t Il3l . Il4 1 1 51 ] is one of the challenging problems 
in heavy quarkonium physics and NRQCD 16]. For 
e + e~ — > J /ipi] c the QCD radiative correction has turned 
out to be essential to g reatly enhance the theoretical pre- 
diction in NRQCD|l7tfl8l]. However, the cross sections of 



other processes, i.e., e H 
by BelleQ 



J/ip{ip(2S))xcO measured 



a[J/ip + Xco] X B*°°[> 2] = (16±5±4) fb, 
^(25)+Xco]xB M [>2] = (17±8±7)fb, (1) 

are also larger than LO NRQCD predictions by about 
an order of magnitude or at least a factor of 5. Here 
B Xa0 [> 2] is the branching fraction for the XcO decay 
into more than 2 charged tracks. Theoretically, two 
studies in NRQCD by Braaten and Lee@] and by Liu, 
He, and Chao[7j showed that, at LO in the QCD cou- 
pling constant a s and the charm quark relative veloc- 
ity v, the cross-section of e + e~ — > J/ilj(ip(2S))xcO at 
= 10.6GeV is only about 2.4 - 6.7(1.0 - 4.4)fb (de- 
pending on the used parameters, e.g., the long-distance 
matrix elements, m c and a s ). So, it is crucial to verify 
that the QCD radiative correction can also greatly en- 
hance <j[J/ijj(ip(2S))xco], before we can claim that the 
large discrepancy between theory and experiment for 
double charmonium production is really resolved. 

However, we do need proof for the validation of fac- 
torization in exclusive production processes involving P- 
wave states at NLO in NRQCD. In fact, nonfactorizablc 
infrared (IR) divergences are found in e.g. the P-wave 
charmonium production in B meson exclusive decays 
such as B — > XcqK[19], in contrast to the factorizable 
S-wave charmonium production in B — > J lib(r\r\K [20|. 
though the IR di verg ence in the P-wave case is m c /mb 
power suppressed 19( . This nonfactorizable feature for 



the P-wave states is a quite general result and is essen- 
tially due to the non-vanishing relative momentum be- 
tween the heavy quark and antiquark in P-wave states. 
So, differing from the double S-wave charmonium pro- 
duction, where factorization can be expected to gen- 
erally hold at NLO in QCD, it is crucial to prove 
the validation of factorization in the special process 
e + e~ — > J/ip(iJj(2S))xc0i where the P-wave charmo- 
nium is involved. Recently, the color transfer has been 
noticed[2lj in associated heavy-quarkonium production, 
e.g., e + e~ — > J/ipcc, where IR divergence appears due 
to soft interactions between the associated c (or c) quark 
and the cc pair of charmonium, and hence breaks down 
factorization at NNLO. Therefore, it is significant to clar- 
ify related problems in QCD for e + e~ — > J/ip(ip(2S))xcO- 
Moreover, NLO QCD corrections are also important in 
understanding heavy quarkonium production at hadron 
colliders (22]]. 

In this paper we will prove the validation of factoriza- 



tion for 



e T e 



J/ipXcO at NLO in QCD, and calculate 



the radiative corrections, while we have already found 
the NLO QCD corrections to e + e~ — > J/^c [l7| and 
e + e~ — * J/ijj + cc[2^] to be large, and increase the cross 
sections by a factor of about 2. All these are at LO in v. 

At LO in a s , J/V'+XcO can be produced at order a 2 a 2 s , 
for which we refer to e.g. Ref [7j. The Feynman diagrams 
are shown in Fig. [T] Momenta for the involved particles 
are assigned as e~ (ki)e + (fc 2 ) — ► J/ijj(2pi) + Xco(2p2)- In 
the calculation, we use FeynArts [24J to generate Feyn- 
man diagrams and amplitudes, FeynCalc [251 ] for the ten- 
sor reduction, and LoopTools [26| for the numerical eval- 
uation of the infrared (IR)-safe one-loop integrals. 

At NLO in a s , there are ultraviolet (UV), IR, and 
Coulomb singularities. We choose renormalization 
schemes the same as in 17], and use D = 4 — 2e dimen- 



sion and relative velocity v to regularize IR and Coulomb 
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FIG. 1: Two of four Born diagrams for e e + — » J/ipXcO- 

singularities. For the box diagrams shown in Fig. [2lBox 
N5, N8, N10 have IR and Coulomb singularities, Box N3 
is IR finite, while the other nine diagrams have IR singu- 
larities. There are some scalar functions that should be 
calculated by hand, which are usually defined as 



Tq[pi ■ ■ .pi-i,m . . .rrn-i] = / -— T 



i(q + pi) 2 - mf 



where p = 0, and T l = C,D,E for i = 3,4,5 respec- 
tively. Some of them were given in (l7j . The others are 
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Here g is the relative momentum of charm quark in Xco, 
r = 16to 2 /s, and Co is the Coulomb and IR divergent 
three point function Co[p c , — Pe, 0, to, to][T3], 
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The IR terms of BoxN5 + N8 + PentagonN10 are canceled 
by counter terms, and the Coulomb singularity is mapped 
into the wave functions. Other IR terms can be separated 
into three point functions Co[p2 + q, — Pi] and Co[p>2 — 
9: — Pi] |22| • And we have 



Co [P2 + q, ~Pi] | q=0 = C [p 2 - q, -Pi] | „ = 



<9C [p 2 + q, ~Pi] 



dq° 



q=0 



dC a [p 2 



q=0 

q, -pi] 



dq c 



(3) 
(4) 



Here Co[l',—l] means Co[l', — I, 0, to, to]. Then we get 
that BoxNl + N4 and BoxN6 + N7 +PentagonN12 are IR 
finite. The IR terms of BoxN9 + N2 + PcntagonNll are 
canceled by vertex diagrams. The UV terms are can- 
celed by counter terms. Then the final NLO result for 



the cross section is UV-, IR-, and Coulomb-finite. De- 
tails of the calculation can be found in a forthcoming 
paper. Since the result is IR finite, e + + e~ — > J/ip + Xco 
is factorizable at NLO in NRQCD factorization. 

The key part for the cancelation of IR divergence in 
our calculation is shown in Fig. [3] The IR term of the 
NLO vertex correction for two charm quarks of momenta 
pi and P2 + q shown in Fig. BKa) is proportional to pi • 
(p2+q)Co[— pi,p2+gi] (see also [2l|,[27[ though the three- 
point function Co was not explicitly written there) ; while 
it is proportional to — pi • (p2 — q)Co[—pi,P2 — q] f° r a 
charm quark and an anti-charm quark of momentum pi 
and P2 — q as shown in Fig. [31(b) . Then for the charm 
quark c(pi) associated with a colorless charm quark pair 
C (P2 + q)c(p2 — q)i the IR term becomes 

MjFtoWpi) + c(p 2 + q)c{P2 - q)] 
«(P2+g) ■ PiC [p 2 +q, -Pi] - (P2~q) -PiC [p2~q, -Pi] 

dC Q p2+g,-pi] I 
dq a 



= 2q a pi a Co [p 2 , -pi] +P2 -Pi ■ 



where it is expanded in powers of the relative momentum 
q at q = 0. (The IR terms between c(p2 + q) and c(p2 — q) 
or between c(pi) and c(pi) are ignored since these cc pairs 
should evolve to bound states at large distances.) 
Following implications can be found from Eq.([5]): 

(1) The IR term is finite when an associated charm 
quark connects with both legs of the S-wave state J/tp, 
where q = can be taken at LO in v, and this corre- 
sponds to the zeroth order i.e. the vanishing O(q ) term 
in Eq.©. 

(2) The IR term becomes divergent when the associ- 
ated charm quark connects with both legs of the P-wave 
state XcOi where the relative momentum q has to be re- 
tained, and this corresponds to the first order i.e. the 
0{q 1 ) term in Eq.©, 



■Pi a s 

S Sir 



1 



~0(q 2 ). (6) 



This nonvanishing IR divergence, which is actually in- 
dependent of the associated quark flavor, is the origin 
for the non-factorizability in many processes involving 
P-wave states, e.g., in B — > XcoK decay 19], and also 
in B decays B — > MiM 2 with Mo being an emitted P- 
wave light meson ( fa, a±, b\...) [28[ when the light meson 
mass e.g. to/ /tob is not ignored in the IR divergent ver- 
tex corrections. In another words, factorization holds up 
only to terms that are m q /mb power suppressed. 

(3) When the associated fermion is an anti-quark c(pi), 
we can get a similar IR divergence by replacing pi with 
— pi in Eq.© and Eq.©. Then by adding together the 
contributions of the associated charm pair c(pi) and c(pi) 
connected with c(p2+q)c(p2—q) , the IR divergence is can- 
celed in the case of P-wave e.g. Xco- Note that since gen- 
erally the associated quark pair c(pi)c(p' 1 ) has pi ^ p[ as 
shown in Fig. [31 the IR cancelation is incomplete and the 
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FIG. 2: Twelve of the twenty-four box and pentagon diagrams 
are for J ftp while two lower ones for \c0- 




FIG. 3: Half of the diagrams for one-loop virtual IR cor- 
rections with two charm quark pairs c(pi)c(p' 1 ) and c(p2 + 
q)c(p2 — q)- The other two diagrams can be obtained by re- 
placing c(pi) with c(p'i). 

divergence still remains. This means that factorization 
in e + e" — > J/ipXcJ can hold at LO in v for the J/tp (i.e. 
Pi = p[) but not hold at NLO in v (pi ^ p[). Based on 
Eq. ([5]) we can draw a general conclusion that the double 
charmonium (including all S,P,D,... wave states) pro- 
duction in e + e~ annihilation is factorizable at NLO in 
a s only on condition that one of the double charmonium 
is an S-wave state in which the quark relative momen- 
tum is ignored. Or, factorization holds up only to terms 
that are to^/s power suppressed. A similar conclusion 
is also obtained recentl y in a more general analysis for 
quarkonium production[29j . 

We now turn to numerical calculations for the cross 
sections of e+e~ — > J/tp(ip(2S))xcO- To be consis- 
tent with the NLO result the values of wave func- 
tions squared at the origin should be extracted from 
the leptonic width of J /ip(ip(2S)) and the two-photon 
width of Xco at NLO in a s (see [l6[ and [I(|), we 
obtain \Ris(0)\ 2 = LOlGeV 3 , |-R 2 s(0)| 2 = 0.639GeV 3 , 
|i?' 1P (0)| 2 = 0.0575GeV 5 . Taking m J/v , = m M2S) = 

m Xc0 = 2m at LO in v, m = 1.5 GeV, = 338MeV, 
we find a s (/j,) =0.259 for (i = 2m, and get the cross sec- 
tions at NLO in a s 

a(e+ + e~ J/i/> + x c0 ) = 17.9fb, 
a(e+ + e- - ^,(2S) + Xc0 ) = 11.3fb, (7) 

which are a factor of 2.8 larger than the LO cross sec- 



3 




Box N4 Box N5 Box N6 




Pentagon N10 Pentagon Nil Pentagon N12 



for e (ki)e + (fe) — * J/V'(2pi)Xco(2p2). Two upper charm legs 



q m = 1.5GeV \ft = 10.6GeV A = 0.338GeV 
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FIG. 4: Cross sections of e + e~ — > J/ip + XcO as functions of 
the renormalization scale fi. 



tions 0.35(4.02) fb for J/ip(ip(2S)). If we use the BLM 
scale [3~0| . we get \iblm = 2.30GeV, a s = 0.291, and 
the corresponding cross sections 8.02(5.08) fb at LO and 
22.8(14.4) fb at NLO. Fig. [4] shows the cross sections at 
LO and NLO as functions of the renormalization scale 
fi, as compared with the Belle and BaBar data. Our 
LO and NLO results compared with experimental and 
other theoretical cross sections are shown in Table I. We 
see the NLO QCD correction enhances the cross sections 
by about a factor of 2.8, despite of existing theoretical 
uncertainties. 

We emphasize again the crucial rule of the associated 
S-wave state J ftp payed in canceling the IR divergen- 
cies in the vertex corrections in e + e~ — > J/ipXcO- At 
LO in v and NLO in a s , the interaction of Xco with the 
charm quark (or antiquark) in the J/ip is individually 
IR divergent, but the sum of that of the charm quark 
and antiquark in the J/ip becomes IR finite. This result 
reflects the fact that the P-wave state Xco behaves as a 
color dipole, which interacts with the color charge carried 
by the charm quark (or antiquark) in the J/ip, but the 
interactions vanish when the charm quark and antiquark 
in the J/ip are combined into a colorless S-wave object 
at LO in v (see Eq.©). The validation of factorization 
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TABLE I: Experimental and theoretical cross sections of 
e + e~ — > J/ip(ifj(2S))XcO at B factories in units of fb. We use 
\Ris{0)\ 2 = l.OlGeV 3 , \R 2S (0)\ 2 = 0.639GeV 3 , \R' 1P (0)\ 2 = 
0.0575GeV 5 , A = 0.338GeV, = 10.6GeV, m c = 1.5GeV, 
and fi = 2m c . The experimental data are the cross sections 
times the branching fraction for XcO decay into more than 2 
charged tracks. But the Belle data of ip(2S) + XcO in Ref. 3] 
correspond to XcO decay into at least 1 charged tracks. 





J/ip + XcO 


V>(2S) + XcO 


Belle a x B x " n [> 2] [1] 


16 ±5 ±4 


17±8±7 


Belle a x B x '"[> 2(0)] [3] 


6.4 ± 1.7 ± 1.0 


12.5 ±3.8 ±3.1 


BaBar a x B Xc0 [> 2] [4] 


10.3 ±2.5^8 




Braaten and Lee [6] 


2.4 


1.0 


Liu, He and Chao [7] 


6.7 


4.4 


Braguta et al. [9] 


14.4 


7.8 


Our LO result 


6.35 


4.02 


Our NLO result 


17.9 


11.3 



at NLO for e + e~ — > J/V'XcO depends crucially on the 
associated S-wave state J/i/>. 

In conclusion, we find that at NLO in a s and LO in v, 
NRQCD factorization holds for the double charmonium 



production 



J/ip(ip(2S))xcO- We get UV and IR 



finite corrections to the cross sections at *Js = 10.6 GeV, 
and the NLO QCD corrections can substantially enhance 
the cross sections with a K factor (the ratio of NLO to 
LO ) of about 2.8; and hence it crucially reduces the 
large discrepancy between theory and experiment. With 
rn = 1.5GeV and /x = 2m, the NLO cross sections are es- 
timated to be 17.9(11.3) fb, which reach the lower bounds 
of experiment. 
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